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We study an electrostatic qubit monitored by a point-contact detector. Projecting an 
entire qubit-detector wave function on the detector eigenstates we determine the precision 
limit for the qubit measurements, allowed by quantum mechanics. We found that this 
quantity is determined by qubit dynamics as well as decoherence, generated by the 
measurement. Our results show how the quantum precision limit can be improved by a 
proper design of a measurement procedure. 
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Rapid experimental progress in monitoring of single quantum systems renewed 
the interest in the quantum mechanical limitation of measurement accuracy 1 . This 
limitation is originated from the uncertainty relation between observables of a mea- 
sured microscopic system, the so-called standard quantum limit of measurement 2 . 
Yet, a single quantum system is not observed directly, but through the interaction 
with a measurement device (detector) . This implies that the quantum limit of mea- 
surement is not related to a measured system only, but rather to an entire system, 
including detector 2 . 

A measurement in the quantum mechanical formalism corresponds to projection 
of the wave function of an entire system on eigenstates of the detector, accessible 
by an observer. This is the so-called projection postulate 3 , analogues to the Bayes 
principle in any probabilistic description. Since the detector and the measured sys- 
tem are interacting, the above projection on the detector states affects the measured 
microscopic system. If the latter is projected on one of its own states, then the mi- 
croscopic system can be measured with any accuracy. However, if this system is 
projected on the superposition of its states, the measurement cannot be precise. Its 
accuracy is given by a size of the corresponding wave packet. 

In this Letter we apply the above described procedure for a determination of 
the precision limit of quantum measurements. As an example we take an electron 
trapped inside a double-dot (electrostatic qubit) and continuous monitored by a 
point-contact detector . The total system can be treated entirely quantum mechan- 
ically, although the detector represents a macroscopic (mesoscopic) device 5 . This 
allows us to obtain the detector eigenstates and then to perform projections of the 
total wave function on these eigenstates without any additional assumptions. 
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Consider a qubit, represented by an electron in the double-dot, placed near 
a point-contact that separated two reservoirs (the point-contact detector), Fig. 1. 
When the first dot, which is far away from to the point-contact, is occupies (Fig. la), 
the current through the point-contact is I\. If the second dot, close to the point 
contact, is occupied (Fig. lb), the current decreases (I 2 < h) due to the elec- 
trostatic repulsion. The entire system can be described by the following tunneling 
Hamiltonian 5 : H = Hpc + Hud + Hint-, where 

Hpc = ^2 E l a \ a l + ^2 E r a t a r + ^2 ^lr( a l a r + H.C.), (1) 
I r l.r 

Hdd = E lC \d + E 2 clc 2 + 0(4ci + c\c 2 ) , (2) 
Hint = ^2 ^irc\c 2 {a\a r + H.c.) , (3) 

l,r 

and SCiir = Q' lr — f2; r . Here a\ r {a>Lr) is the creation (annihilation) operator for 
an electron at the level / or r in the left or right reservoir, and c\ 2 {c\_ 2 ) is the 
same operator for the electron inside the double-dot. fli r is the hopping amplitude 
between the states I and r of the reservoirs, and f2 is the hopping amplitude between 
the states E\ and E 2 of the qubit. For simplicity we consider electrons as spin- less 
fermions. The interaction term Hint generates variation of the hopping amplitude, 
(50; r = fl'ir — f2; r , resulting in a decrease of the detector current from I\ to I 2 , 
Fig. 1. 
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Fig. 1: The electrostatic qubit monitored by the point-contact detector. The de- 
tector current decreases when the left dot of the qubit is occupied. Here n denotes 
the number of electrons which have arrived at the right reservoir by time t. 



The wave function describing the entire system can be written as 



!*(*)> = 



h(t)c\ +Y2bu r (t)c\alai + ^ &l(i'rr'(*) c l4 a r' a i a i' 
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|0), (4) 



where b a (t) is the amplitude of finding the system in the state a determined by 
the corresponding creation and annihilation operators. These operators act on the 
initial ("vacuum") state, |0). For simplicity we assume that the reservoirs are ini- 
tially at zero temperature and filled with electrons up to the Fermi energies ^l,r, 
respectively. All the amplitudes b a (t) can be obtained from the time-dependent 
Schrodinger equation, d t \^(t)) = H\^(t)}. 

Now we project the total wave function (4) on the detector eigenstates. There 
exists an uncertainty however, about a choice of the eigenstates, since different de- 
tector variables can be recorded 6 . In fact, the point-contact detector is not recorded 
directly, but via another readout device ("pointer"). The latter can single out a 
particular set of the eigenstates through a coupling with the corresponding detec- 
tor variables. Yet, in this work we do not extend our system by including such a 
pointer in the Schrodinger equation. We assume instead that the detector states 
are directly accessible to an "observer" . Let us examine different alternatives of the 
measurement. 

Consider first the measurement of number of electrons (n) in the right reser- 
voir (the accumulated charge), Fig. 1. This implies that the wave function, |\Jj(t)), 
Eq. (4), is projected on the eigenstates of the operator N = J2 r a t a r, which can be 
written as 



\n) = a. 



|0>. 



(5) 



(Note, that the state \n) is strongly degenerate: \n) = |n, a), where a corre- 
sponds to a particular configuration of n electrons in the collector). Thus, | — ► 
J\f- 1 / 2 P n \^(t)}, where P n = J2 a \n,a){n,a\ is a projection operator and TV is a 
normalization factor. One finds from Eqs. (4) and (5) that 



P n |*(t)> = £ 




10} (6) 



where the corresponding normalization factor TV = J2i r [\bu...r{t)\ 2 + \b2i... r (t)\ 2 ] 
is a probability of finding n electrons in the collector by time t. Eq. (6) shows that the 
measurement leaves the qubit in a linear superposition of its two states. Therefore 
one cannot determine the state of the qubit by measuring the number of electrons 
in the right reservoirs. 

Consider now the measurement of electric current in the right reservoir. The 
latter is given by a commutator of N with the total Hamiltonian of the system, 
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I = i[H, N] (we choose e = 1). Using Eqs. (l)-(3) we obtain 

I = iy^(Q; r + 5n lr c\c 2 ){a\a r - a\.a{) = ^ f ir (7) 

Z,r l,r 

The eigenstates of the energy resolved-current operator, Ii r \Ij^(q)) = I^(q)\Ii^(q)) 
are 

|J±(g)) = -^(at Ol ±i)ct|0), (8) 

where g = 1,2 denotes the qubit state, c||0) and /^(<z) = ±[flj r + (5 — l)£fi; r ]. 
Respectively, the eigenstates of the total current, |7), are given by a product of 

It follows from Eqs. (7), (8) that by measuring the energy-resolved current Ii r 
(or the total current, I), one projects the wave function (4) on a certain state of 
the qubit, corresponding to an observed value of the current. This implies that the 
qubit position can be determined with any accuracy, in principle, by monitoring 
directly the detector current (via its magnetic field 7 ). 

If, however, such a direct measurement of the detector current cannot be per- 
formed, one can determine it indirectly, via a variation of the collector charge. Let 
us assume that we recorded no electrons in the collector at time t. As a result, the 
entire system is projected to the state \no): \ip(t)) — ► P no \ip(t)), Eq. (6), which is 
an eigenstate of the operator N, Eq. (5). (We omitted the index of degeneracy a). 
Next we detect the accumulated charge n at the time t + At. The final state of the 
detector (up to the normalization factor) is 

P n e-* HAt P no \y(t)) = \n){ Vn (At)\n ){n \* (*)) , (9) 

where \ip n (At)) = cxp(iHAt)\n) is an eigenstate of the operator N(At) — 
cxp(iHAt)Nexp(—iHAt), corresponding to an eigenvalue n. This operator can be 
expanded in powers of At 

N(At) = N + IAt + i[H, 1]^—^— + • • • , (10) 

where the current I is given by Eq. (7). Thus the time-dependent operator N(At) 
includes the qubit position operator, cjc 9 , in contrast with N = N(0). As a result 
the projection on the eigehstates of N(At), Eq. (9), could determine the qubit 
position. 

Let us take small At ("measurement time") in Eq. (10) such that N(At) ~ 
N + IAt. If and I commute, then the eigenstates of N(At) would be a product 
of eigenstates of these operators: \ip n (At)) — \n',I(q)), where q = {1,2} denotes 
the qubit state and n = n! + I(q)At. It follows from Eq. (9) that n' = no and 
I(q) = iAn(q) = An/ At, where An = n — no. As a result the qubit is projected in 
the state q corresponding to the variation of the collector charge, An. 

In fact, the operators N and I do not commute. In this case it is only the average 
current, I& n = An /At which can be determined from ensemble measurements of 
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An, where its dispersion, 5I\ n = [(An) 2 — (An) 2 ] 1 ' 2 /At, diverges as (At) -1 / 2 for 
At — > 0. The latter restricts the accuracy of the qubit measurements, respectively. 

The measurement accuracy, however, increases by increasing the measurement 
time, since 51& n {q) — ► for At — > oo 8 . Yet, this can be done only if the qubit is 
not "moving", (i.e. the hopping amplitude £1 = 0, Fig. 1). Then [H, I] = 0, so that 
the higher order terms in the expansion (10) vanish. If it is not the case (O ^ 0), 
the current it driven by the qubit, so that [H, I] and the higher order commutators 
in Eq. (10) are not zero. The average contribution from these terms to I\ n , which 
we denote as 8iI\ n (/S.t), increase with At. This suggests that the quantum limit of 
the qubit measurement is determined by the optimal measurement time (At) which 
minimizes the total error, 

[<5 2 /A„(At)] 2 = [SI An (At)} 2 + [SilAniAt)} 2 . (11) 

In order to perform this procedure we introduce the density matrix cr^™, (t) = 
(n, q'\^(t))(^>(t)\n, q), where the wave function |^(t)) is given by Eq. (4). It was 
demonstrated in 5 ' 9 that for large large bias voltage, V — u £ — p R (Fig. 1), the 
time-dependent Schrodinger equation, dt\^(t)) = H\&(i)}, can be reduced to the 
following Bloch-type rate equations for the density matrix a^™} (t) by assuming weak 
energy dependence of the transition amplitudes (Q; r = Cl, Q' lr — Cl') 5 , 

= -D x a<$ + D^ 1 ^ + iQ{a[f - 4°) (12) 
( " } - -D 2 *%> + D 2 *&-V in(a[? ,W) (13) 



, 22 



■IT? r? \ (") i •O/' (") ( n )\ D\ + D 2 (n) /„ „ 0/2 (n-l) /1 .n 

'12 - A E 2 - Ei)<j\ 2 > + iQ.{(t\ 1 ' - o\ 2 >) o\ 2 > + (D 1 D 2 ) 1 a\ 2 {14) 

where D\ 2 = T\^Y and 7i ; 2 is the transmission probability of the barrier: T\ = 
(2tt) 2 0, 2 plPr and T 2 = (2tt) 2 (CI') 2 plPr, where pl,r is the density of states in the 
left (right) reservoir, Fig. 1. 

Solving Eqs. (12)-(14) one can find all quantities needed for evaluation of 5 1 An 
and 8\I\ n - For instance, in order to evaluate the average value of An and its dis- 
persion we have to solve these equations with the initial condition n = no- It follows 
from Eqs. (12)-(14) that the density matrix cr^(t) depends only on An = n — n . 
Thus we can take n = and An = n. Then the average values An = n and 
(An) 2 = n 2 are given by 

n{t) = ^nP n {t)=n 11 {t)+n 22 {t), (15) 

n 

^(t) = ^n 2 P„(t)=<(t)+^(t), (16) 

n 

where P n (t) = {t)+&22 (*) i s the probability of finding n electron in the collector. 
Multiply Eqs. (12)-(14) by n and sum over n one finds 

nn = -D1011 + iCl(ni2 — mi) (17) 
«22 = D 2 a 22 - iH(ni2 - n 2 i) (18) 
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ni2 = i{E 2 - Ei)ni2 + iCl(nu - n 22 ) - "y^i2 + (DiD 2 ) 1/2 <t 12 , (19) 
where <J qq >(t) = J^n^qq'W i s the qubit density matrix, and Td = (\/Di — is 



99 

the decoherence rate 5 . 

Similarly multiplying Eqs. (12)-(14) by n 2 and sum over n one obtains 

nf^ = 2Dinu + DiCTii + ifi(nf^ - nfj (20) 
n2 2 = 2D 2 n 22 + D 2 a 22 - i£l(n\ 2 - ra^i) (21) 
n£ = - -Ei)n£ + - "ra) ~ y ™?2~ + (^i^ 2 ) 1/2 (2ni 2 + <t 12 ) . (22) 

The qubit density matrix <J qq >(t) can be easily found from Eqs. (12)-(14) by 
tracing it over n, 

an = ifl(cri2 - CT21) (23) 

0"12 = «e2lCTl2 + ifi(2(Tll - 1) - y(Tl2 , (24) 

and a 22 (t) = 1 — (Tn(t). This quantity a qq i(t), determines the detector average 
current, 7(t) = (*(i)|/|*(t)). Indeed, it follows from Eqs. (4), (7) that 

7(t) = DiTiuit) + D 2 ji 22 (t) = D 2 + AD an (t) , (25) 

where AD = D\ — D 2 is an average "signal". Obviously, I (At) =n/At = I n (At) 
for small At. The dispersion of I n (At) can be found from Eqs. (16), (20)-(22), 

(5I n f = ^(At)l(At) 2 -t n ~ I n (Q)/At . (26) 

As expected, SI n (At) diverges as 1/VAt for At — » 0. 

Let us evaluate the contribution from higher order terms in the expansion (10), 
which generate variation of the average current, 5\I n . (In our case this variation is 
produced by the qubit only). Therefore, one can write 8\I n (At) = \I n (At) — I n (0)\. 
Using Eq. (25) we evaluate this quantity as 

5il n = AD |<7ii(Af) — crii (0)| , (27) 

where <rn(t) is obtained from Eqs. (23)-(24). For instance, for aligned levels, Ei — 
E 2 , and ctu(0) = 1, 



Sil n = — 



c 3~ [cos(wAt) + 77 sin(wAt)] - 1 , (28) 



where r) = Td/Auo and uj = 20-^/1 — (Td/8£l) 2 is the Rabi frequency. As expected, 
Sil n -» for n -> 0. 

Eqs. (26) and (28) allow us to evaluate the optimal At by minimizing [<5 2 /„(At)] 2 , 
Eq. (11). We take for simplicity AD < D, where D = (Di + D 2 )/2. As a result 
T d ~ (AL>) 2 /4D. We first consider weak distortion of the qubit, r d /8 < Cl. Then 
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expanding Eq. (28) in powers of At we easily obtain for the optimal measurement 
time and for the corresponding measurement limit 



In order to observe Rabi oscillations of the qubit in a single measurement run one 
needs <C AD, at least. It follows from Eq. (29) that this condition corresponds 
to £1 <C 2Td- This, however cannot be combined with the condition of weak qubit 
distortion, used in Eq. (29). Hence, one cannot observe Rabi oscillations in a single 
run, but only in an ensemble average of different runs. 

Consider now large decoherence limit, 1^/8 » f2. Then the qubit is strongly af- 
fected by the detector. As a result, the electron stays in the same dot for a long time, 
t ~ Td/8fl 2 ("quantum Zeno" effect). Indeed, one finds from Eqs. (23)-(24) that 
(Til = [1 + cxp(— 8V, 2 t/Td)]/2 for t 3> Respectively, the optimal measurement 

time and the measurement limit are given by 



In contrast with the previous case, Eq. (29), the measurement time At increases 
with Td- This is not surprising since large decoherence generated by the detector, 
localizes the qubit for a long time. Therefore it behaves as a static one so that the 
measurement time increases. 

It follows from our arguments that the quantum precision limit depends on a 
particular set of the detector observables which the total wave function is projected 
on. In this Letter we discussed two alternative sets related to charge and current 
states of the point-contact detector. It was demonstrated that single projection on 
charge states cannot measure the qubit state. However, two consecutive projections 
of the entire system on the charge states can measure the qubit state, but only with 
a limited accuracy. On the other hand, direct projection on the current state of the 
detector can determine the qubit state with absolute precision, in principle. If this 
could be realized, one would arrive to the Zeno paradox 10 , i.e. to complete freezing 
of a system as a result of continuous measurement. 

This shows a necessity of including a "pointer" in the total Hamiltonian coupled 
with current states of the detector (von Neumann hierarchy 3 ). In this case two 
consecutive projections of the total wave function on the pointer states, used for a 
determination of the detector current, would restrict the measurement accuracy in 
a total analogy with the previous case. If the pointer is coupled with charge states 
of the detector, it obviously cannot decrease the quantum measurement limit found 
in our calculations. We assume also that the pointer cannot essentially increase 
this limit, since otherwise the von Neumann hierarchy of measurements would not 
converge. This problem, however needs a special attention. 

Our final results on quantum limit of measurement involved only average quan- 
tities, which were obtained without any explicit resort to the projection postu- 




(29) 




(30) 
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late. Indeed, the latter is related only to a single measurement. Nevertheless, as 
we demonstrated in this Letter, the use of the projection postulate was very use- 
ful in a determination of quantum limit of measurement. In particular, it clearly 
displayed the detector variables which would allow us to measure a microscopic 
system with maximal accuracy. Such variables are usually represent commutators 
(time derivatives) of operators describing the detector states. 

An appropriate choice of this variable depends of a particular measurement 
apparatus. For instance, for a single electron transistor (SET) detector, one needs to 
use the second commutator ("acceleration") of the accumulated charge. In contrast 
with the point-contact detector, the projection to current states of the SET would 
not determine the qubit state. The measurement of the charge "acceleration" can be 
very useful if the corresponding operator would commute with the charge operator. 
In this case one can design an appropriate procedure of projecting on the charge 
states at different times which would diminish the quantum measurement limit. 
This however must be a subject of a separate investigation. 
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